X 


«I>-A053  567 


UNCLASSIFIED 


MISSISSIPPI  UNIV  UNIVERSITY  DEPT  OF  ELECTRICAL  EN6IN--ETC  F/6  20/3 
A STUDY  OF  CHARGE  PROBE  RESPONSE  ON  A JUNCTION  OF  UNEQUAL  RADII—ETCCU) 
SEP  76  S SINGARAYAR.  B H DUFF  AFOSR-75-2862 


AFOSR-TR-78-0772 


NL 


U Ht  i 


EPORT  DOCUMENTATION  PAG]; 


2.  GOVT  *.CCrSSICN 


^ ^DY  OF  SPARGE  JROBE  ^SPONSE  ON  A JUNCTION 
tfif^EQUAL  SlDII  mN  WI^S^  "S- 


Kf/TD  INSTKiJCTIONS 
EFGKC  COMPl  rriNff-KOHM 


-IPTKT'S  CAT  AwOJ'  NUMBER 


9.  PERrORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Electrical  Engineering 
University  of  Mississippi 
University,  Mississippi  38677 


It.  CONTROLLING  OFFICE  NAME  AND  AIVBiBI 

AFOSR/NF  f\\j 

Bolling  AFB,  Bldg. #410 llj/j 
Wash  DC  20332  ' 


14.  MONITORING  AGENCY  NAME  & ADORES 


2-e.r^  IL 


15.  SECURITY  CLASS.  ( 


Unclassified  


15«.  OECLASSIFICATION/OOWNCR  AGING 
SCHEDULE 


16.  distribution  statement  ro/ n^porO 


Approved  for  public  release; 
distrlbutloa  unliffllted* 


17.  distribution  statement  (ot  lh«  mbttfct  tnt^fd  in  Block  1/  ditl^rent  from  Report) 


16.  supplementary  notes 


19.  KEY  WORDS  ^Conllnuo  on  revetee  eide  If  neceesary  and  iden(//y  by  block  number) 

Charge  Probe 

Response  on  Different  Radii  Cylinders 
Experimental  Study 
Theoretical  Study 


25r''*'A^^TRACT  (Continue  on  revoraa  aide  If  neeeaaary  and  idonfllx  by  block  number) 

The  experimental  measurement  of  charge  density  induced  on  conducting  surface 
often  makes  use  of  a short  monopole  as  the  charge  probe.  Measurements  on  struc- 
tures containing  junctions  of  unequal  radii  cylinders  have  been  hampered  by 
incomplete  knowledge  of  the  probe  characteristics.  An  investigation  of  charge 
probe  response  as  related  to  the  radius  of  the  cylinder  on  which  it  is  mounted 
is  reported  here.  An  experimental  determination  of  the  ratio  of  probe  response 
when  mounted  on  two  different  radii  cylinders  has  been  conducted  using  a coaxial 


DD  /.Tun  1473 


EDITION  OF  I NOV  6r  IS  OBSOLETE 


UncLnssi  i 3 

security  CLASSIFIC*’  I or  T .IS  PAGE  flUi.n  Dmlm  Enttfd) 


f \ 

Unclassified 

SECur)IT>  classification  of  this  PAGEfltTicn  Dmt»  F.nlarod) 

' line  geometry  with  a step  change  in  radius  of  the  inner  conductor.  In  addition, 
a theoretical  solution  in  numerical  computations  of  the  probe  response  in  rela- 
tion to  the  radius  of  the  cylinder  upon  which  it  is  mounted  has  been  conducted. 
The  results  of  these  studies  are  presented  in  the  following  report. 


I 

j 


A STUDY  OF  CHARGE  PROBE  RESPONSE  ON  A JUNCTION 
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Abstract 


The  experimental  measurement  of  charge  density  induced  on  conducting 
surfaces  often  makes  use  of  a short  monopole  as  the  charge  probe.  Measure- 
ments on  structures  containing  junctions  of  unequal  radii  cylinders  have 
been  hampered  by  incomplete  knowledge  of  the  probe  characteristics.  An 
investigation  of  charge  probe  response  as  related  to  the  radius  of  the 
cylinder  on  which  it  is  mounted  is  reported  here.  An  experimental  deter- 
mination of  the  ratio  of  probe  response  when  mounted  on  two  different 
radii  cylinders  has  been  conducted  using  a coaxial  line  geometry  with  a 
step  change  in  radius  of  the  inner  conductor.  In  addition,  a theoretical 
solution  in  numerical  computations  of  the  probe  response  in  relation  to 
the  radius  of  the  cylinder  upon  which  it  is  mounted  has  been  conducted. 

The  results  of  these  studies  are  presented  in  the  following  report. 
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Chapter  I 
INTRODUCTION 

The  measurement  of  the  current  and  charge  distributions  on 
antennas  and  scatterers  is  important  in  order  to  gain  an  under- 
standing of  their  electromagnetic  behavior  and  to  aid  in  developing 
and  checking  of  theoretical  solutions.  The  basic  techniques  for 
making  these  measurements  have  been  established  for  some  time  and 
used  by  many  investigators.  However,  during  a recent  experimental 
project,  to  be  described  below  it  became  obvious  that  more  accurate 
knowledge  of  the  response  characteristics  of  charge  probes  was 
necessary.  The  resulting  investigation  of  charge  probe  response 
forms  the  basis  of  this  thesis. 

The  experimental  project  which  leads  to  the  investigation  re- 
ported here,  concerned  the  measurement  of  charge  distribution 
Induced  on  structures  containing  junction  of  thin  wires  of  unequal 
radii.  The  purpose  of  that  investigation  was  to  experimentally 
determine  the  appropriate  condition  on  charge  which  should  be 
applied  in  theoretical  solution  of  this  type  of  structure. 

A theoretical  analysis  and  discussion  of  the  junction  has  been 
presented  by  T.T.  Wu  and  R.W.P.  King  [3].  The  probe  used  for  the 
measurement  of  charge  distribution  consisted  of  a very  short  mono- 
pole mounted  on  the  surface  of  the  cylindrical  structure.  For 
the  experimental  determination  of  the  charge  condition  at  the 
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junction,  an  accurate  knowledge  of  the  probe  response  as  a function 
of  the  radius  of  the  cylinder  upon  which  it  is  mounted  is  necessary. 
A charge  probe  mounted  on  a cylindrical  structure  has  been  analysed 
theoretically  by  R.W.P.  JCing  [1]  and  experimentally  by  Whiteside  [2]. 
However,  these  investigations  do  not  provide  sufficiently  accurate 
solutions  for  the  probe  response  as  a function  of  the  radius  of  the 
cylindrical  structures.  An  independent  experimental  and  theoretical 
investigation  of  the  charge  probe  response  has  been  conducted,  the 
results  of  which  are  presented  in  this  thesis. 

Chapter  II  deals  with  the  measurement  technique  by  which  the 
probe  calibration  factor  can  be  determined.  This  probe  calibration 
factor  has  been  determined  for  3 different  probe  heights  for  2 
Junctions  of  thin  wires  of  5/16"  to  1/8"  radii  and  3/16"  to  1/8" 
radii.  To  make  sure  that  the  probe  calibration  factor  is  the  same 
for  all  charge  conditions  at  the  junction,  the  probe  calibration 
factor  has  been  determined  for  3 different  charge  conditions  namely 
maximum,  minimum,  and  a condition  in  between  maximum  and  minimum. 

The  experimental  data  was  analysed  by  using  modified  Prony's  method. 

A numerical  solution  for  the  current  induced  in  the  probe 
has  been  carried  out  using  moment  method  techniques  [4,5].  The 
analysis  models  the  probe  as  a monopole  receiving  antenna  perpen- 
dicular to  an  infinite  plane  conductor,  but  Includes  the  incident 
E field  variation  which  is  to  be  expected  when  the  probe  is  mounted 
on  a cylinder.  Chapter  III  deals  with  the  methods  used  for  this 
numerical  analysis. 
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Finally,  in  Chapter  IV,  the  experimental  data,  the  results  and 
the  conclusion  are  presented. 


CHAPTER  II 


MEASUREMENT  TECHNIQUE 

(1)  Introduction 

The  measurement  of  the  surface  charge  density  on  a conductor  may 
be  related  to  a measurement  of  E field  by  using  the  boundary  condition, 

n.epE  = n (2.1) 

where  is  the  relative  permittivity,  n is  the  unit  normal  vector  to 
the  surface  of  the  test  structure  under  consideration  and  n is  the 
surface  charge  density.  It  is  therefore  possible  to  obtain  the  surface 
charge  density  by  using  a sensor  whose  output  is  proportional  to  the 
, normal  component  of  the  E field  at  the  surface.  The  most  common  method 

uses  a short  monopole  receiving  antenna,  as  shown  in  Fig.  (2.1), which 

is  perpendicular  to  the  surface  of  the  test  structure.  The  signal 
[ induced  in  the  load  is  however  related  to  the  E field  along  the  entire 

1 length  of  this  monopole  rather  than  simply  the  field  at  the  surface. 

! If  the  normal  component  of  the  E field  is  constant  over  the  length  of 

! the  probe  there  is  no  problem.  But  generally  this  is  not  the  case. 

In  order  to  relate  the  monopole  probe  response  to  the  charge  density 
on  the  test  structure  it  is  therefore  necessary  to  have  some  knowledge 
of  the  manner  in  which  E varies  over  a short  distance  away  from  the 
surface. 
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The  present  work  considers  only  test  structures  composed  of 

electrically  thin  conducting  cylinders.  The  field  produced  by  a 

uniformly  charged  infinite  cylinder  is  known  to  vary  as  If  the 

test  structure  (a  cylindrical  antenna)  is  of  electrically  thin  radius 

and  the  probe  is  electrically  short,  then  E over  the  length  of  the 

probe  is  controlled  primarily  by  the  local  charge  distribution.  The 

field  then  becomes  essentially  the  same  as  that  of  an  infinite  cylinder 

and  E is  very  nearly  proportional  to  — . This  assumption  is  not  valid 

P 

near  discontinuities  of  the  conducting  surface,  however.  The  effect 
of  this  variation  on  the  response  of  the  probe  depends  on  the  radius 
of  the  test  structure.  To  illustrate  this  effect.  Fig.  (2.2)  presents 
plotted  as  a function  of  (p-a)  for  three  different  cylinder  radii. 

For  a given  probe  height,  the  variation  of  E^  over  the  length  of  the 
probe  is  less  when  the  probe  is  mounted  on  a large  radius  cylinder. 


The  effect  of  the  — variation  may  be  reduced  by  reducing  the 
probe  height.  However,  the  signal  induced  across  the  load  is  also 
reduced  and  so  the  signal  to  noise  ratio.  A second  effect  of  changing 
the  radius  of  the  test  structure  is  that  the  image  plane  for  the 
monopole  probe  is  of  different  radius  of  curvature.  The  fact  that  the 
charge  probe  is  of  finite  size  and  it  responds  to  fields  away  from  the 
surface  of  the  test  structure  makes  it  necessary  to  determine  a probe 
calibration  factor  as  a function  of  cylinder  radius  in  an  exact  manner. 

To  determine  the  probe  calibration  factor  experimentally,  a system 
containing  cylinders  of  different  radii  and  for  which  the  variation  of 


the  charge  density  as  a function  of  cylinder  radius  can  be  predicted 
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is  needed.  The  system  chosen  was  a coaxial  line  with  step  change  in 
the  radius  of  the  inner  conductor.  This  system  Is  shown  In  Fig.  (2.3) 
and  has  been  analysed  by  Marcuvltz [6] . The  TEM  mode  E field  dis- 
tribution of  the  above  system  is  shown  in  Fig.  (2.4).  The  E field  is 
uniform  and  proportional  to  inside  the  coaxial  line  except  near  the 
junction.  The  local  effect  near  the  junction  can  be  replaced  by  a 
shunt  capacitance  .so  that  the  Kirchoff's  current  law  is  satisfied. 

At  the  junction  the  voltage  is  continuous  making  it  easy  to  analyse 
the  charge  condition  at  the  step  in  the  following  way. 


(li)  Analysis  of  Coaxial  Transmission  Line 

Fig.  (2.5)  illustrates  a coaxial  transmission  line  with  a TEM 
mode  propagation.  In  cylindrical  coordinates  the  two  dimensional 
Laplace  equation  is. 


1 

P 


3 

8p 


P 


= 0 


(2.2) 


For  a potential  function  independent  of  angular  coordinate  (f), 
eq.  (2.2)  reduces  to 


1_ 

P 


3p  3p^ 


0 


(2.3) 


The  solution  of  equation  (2.3)  may  be  expressed  as, 

= Cj^  i-n(p)  + c^  (2.4) 

Imposing  the  boundary  conditions  at  p = a and  $ = 0 at  p = b, 

eq.  (2.3)  yields  the  final  solution  for  the  potential  as. 
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♦ - Vq  ln(  ^ )/£n(  ) 

E field  of  TEM  mode  in  +z  direction, 

E = - 4>  e"^*“ 


where  k is  the  propagation  constant. 
Using  eq.  (2.5)  E may  be  expressed  as. 


E 


iln(^) 


/V 

a .t 
_P  g-jkz 

p 


(2.5) 


(2.6) 


(2.7) 


The  surface  density  of  change  on  an  inner  conductor  of  radius  a^^, 
at  z=0  may  be  obtained  from  eq.  (2.7)  and  the  boundary  condition 
eq.  (2.1)  as. 


^1  = 


a,£n(^)  0 


(2.8) 


Equation  (2.8)  can  be  written  as. 


where  is  the  capacitance  per  unit  length  and  q^^  = 2TTa^r|j^  is  the 
linear  charge  density.  The  same  relationship  applies  to  the  inner 
conductor  of  radius  a2- 

It  must  be  remembered  that  the  above  analysis  assumes  that  only 
TEM  mode  fields  exist  and  that  the  junction  is  represented  by  a lumped 
shunt  capacitance.  The  effect  of  the  step  change  of  radius  is  actually 
distributed  over  some  short  distance  on  either  side  of  the  junction  with 
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an  associated  charge  density  which  differs  from  that  of  a TEM  mode  as 
indicated  in  Fig.  (2.4)(b). 

The  voltage  is  continuous  at  the  junction  but  because  of  the 
different  radii,  a difference  in  the  magnitude  of  charge  density  exists 
at  the  junction.  The  ratios  of  linear  and  surface  charge  densities  at 
the  junction  are  given  by 


An  (^) 

^ 

An  (^) 
®1 


(2.9)(a) 


An  (— ) 

h ^2 

An  A 
®1 


(2.9)(b) 


(ill)  Description  of  the  Measurement  System 

A block  diagram  of  the  measurement  system  is  given  in  Fig.  (2.6). 
The  commercial  Instruments  used  in  the  system  are  given  in  the  block 
diagram  Itself. 

The  power  at  300  MHz  is  fed  to  the  amplifier.  The  amplifier  sig- 
nal is  fed  to  the  coaxial  system  through  a balanced  feed  to  avoid 
higher  order  mode  propagation.  Matching  of  the  amplifier  output  to 
the  coaxial  system  is  done  through  a single  stub  tuner.  The  reference 
signal  to  the  vector  voltmeter  reference  channel  is  taken  from  the 
amplifier  output  with  proper  attenuation.  The  signal  to  the  frequency 
meter  is  also  is  taken  from  the  amplifier  output  with  proper  attenua- 
tion. The  signal  from  the  charge  probe  is  fed  to  the  test  channel 


of  the  vector  voltmeter.  The  power  level  is  monitored  by  the  broad- 
band sampling  meter  to  compensate  for  the  drift  In  power  level. 

The  magnitude  or  the  phase  of  the  signal  Is  switched  to  the 
digital  voltmeter  by  means  of  the  relay  actuator.  The  calculator 
samples  this  output  5 to  10  times  and  takes  the  average  to  minimize 
the  effect  of  the  minute  fluctuations  in  the  probe  signal.  The  magni- 
tude and  phase  quantities  are  stored  In  the  calculator  for  later 
plotting  and  for  further  use. 

The  stepped  radius  coaxial  line  was  constructed  from  cylindrical 
brass  tubing.  The  three  different  inner  conductor  radii  0.3125" 
(0.79375  cm),  0.1875"(0. 47625  cm),  and  0.125" (0.3175  cm)  were  selected 
as  being  available  readily.  Two  inner  conductor  systems  were  con- 
structed one  composed  of  0.3125"  and  0.125"  radius  tubing,  the  other 
of  0.1875"  and  0.125"  radius  tubing.  A tube  of  1.0"  inside  diameter 
was  used  as  the  outer  conductor.  The  space  between  Inner  and  outer 
conductor  was  filled  with  styrofoam.  A slot  of  1/19" (0.13368  cm)  was 
cut  axially  in  the  inner  conductor  tubing  to  allow  movement  of  the 
charge  probe  along  the  structure.  The  probe  carriage  was  designed 
to  slide  inside  of  the  inner  conductor  tube  and  extend  through  the 
slot  to  the  surface  of  the  tube  so  that  the  structure  is  electrically 
continuous.  The  charge  probe  was  constructed  on  a probe  carriage 
using  a semi  rigid  coaxial  cable.  The  charge  measurement  was  con- 
ducted on  two  different  radii  tubes.  If  two  different  probes  had  been 
used  it  would  have  been  necessary  that  they  be  constructed  Identically. 
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Otherwise,  it  would  have  become  necessary  to  calibrate  them  indepen- 
dently. To  avoid  this  problem,  a single  charge  probe  with  two  different 
probe  carriages  was  used.  The  wall  thickness  of  the  tubes  had  to  be 
taken  into  consideration  while  making  the  probe  carriage.  For  example 
for  the  0.3125"  radius  tube  wall  thickness  was  0.06" (0.1524  cm)  and  the 
0.125"  radius  tube  was  0.04".  The  probe  carriage  fit  exactly  inside 
of  the  tube  and  through  the  slot  to  the  surface  so  that  the  charge 
probe  was  effectively  mounted  on  the  surface  of  the  inner  conductor 
tube.  The  cable  used  for  making  the  probe  was  UT  20  manufactured  by 
Uniform  Tubes,  Inc.,  with  the  following  specifications. 

Impedance:  50  Ohms 

Outside  diameter:  0.023" (0.05842  cm) 

Center  conductor  diameter:  0.0045" (0.01143  cm) 

The  outer  conductor  of  the  cable  was  soldered  to  the  probe  carriage  and 
cut  at  the  surface.  The  dielectric  insulator  on  the  inner  conductor 
was  removed  so  that  the  inner  conductor  extended  beyond  the  surface  to 
form  the  monopole  antenna.  It  was  necessary  that  this  monopole  be 
straight  and  perpendicular  to  the  test  structure.  The  signal  from  the 
probe  was  carried  by  a coaxial  cable  contained  inside  of  the  inner 
conductor  tube.  Care  was  taken  to  see  that  stray  pickup  did  not  affect 
the  signal  by  shielding  the  signal  cables.  A 1/8"  diameter  tube 
attached  to  the  probe  carriage  was  used  both  to  protect  the  signal 
cable  and  to  provide  a means  of  positioning  the  probe. 


The  charge  density  measurement  was  made  at  intervals  of  1 cm.  The 
actual  E field  distribution  near  the  step  due  to  local  effect  made  the 
measured  result  unreliable  in  this  region.  Moreover  the  charge  probe 
was  mounted  at  the  center  of  the  probe  carriage  which  made  it  impos- 
sible to  reach  the  step  in  the  large  radius  tube.  For  analysis  of  the 
data,  a few  data  points  were  eliminated  near  the  junction  to  eliminate 
the  local  effect  and  the  data  was  analysed  by  using  a modofied  Prony's 
method  to  fit  curves  to  the  TEM  mode  distribution.  This  procedure 
will  be  discussed  separately. 

Measurements  were  made  for  probe  heights  of  0.10" (0.2540  cm)  to 
0.06"(0.1524  cm)  at  an  interval  of  0.02"{0.0508  cm).  The  charge 
condition  at  the  junction  was  varied  by  using  the  movable  short  and 
the  experiment  was  carried  out  for  three  charge  conditions  at  the 
junction,  namely  charge  maximum,  charge  minimum  and  a case  in  between 
maximum  and  minimum.  The  experiment  was  repeated  for  consistency. 

Even  though  frequency  drift  was  observed,  the  change  in  electric  length 
of  the  test  structure  was  of  only  of  the  order  of  + O.OOIA  which  was 
considered  negligible.  Any  change  in  the  power  output  of  the  oscil- 
lator was  compensated  for  by  monitoring  the  oscillator  output  using  a 
broadband  sampling  meter. 

Care  was  taken  while  making  the  probe  to  see  that  the  top  of  the 
probe  carriage  fits  flush  with  the  test  structure.  To  achieve  good 
electrical  contact  between  the  probe  carriage  and  the  test  structure, 
silver  conducting  grease  (Eccoshield,  manufactured  by  Emerson  & Cuming) 
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was  applied  along  the  sides  of  the  probe  carriage  and  the  sidewalls 
of  the  probe.  The  signal  carrying  probe  tube  was  coated  with  an 
absorber  (Eccosorb  by  Emerson  & Cuming)  to  avoid  the  possibility  that 
the  probe  tube  would  become  the  center  conductor  of  a TEM  structure 
Inside  of  the  Inner  conductor  of  the  test  structure  and  give  rise  to 
propagation  and  unwanted  resonances. 

After  obtaining  the  extrapolated  value  of  the  probe  reading  at  the 
junction  the  probe  calibration  factor  was  obtained  as  described  in  the 
following  section. 

(iv)  Determination  of  Probe  Calibration  Factor 

Let  and  P2  be  the  probe  reading  obtained  by  extrapolating  the 
measured  data  to  the  junction  of  the  large  and  small  tubes  of  radii  a^^ 
and  Let  and  N2  be  the  normalizing  factors. 

Pj  = (2.10) 


^2  * ^2^2 


(2.11) 


therefore 


!i!2 

P2 


from  eq.  (2.9)(b)  the  ratio  of  — is  given  by, 


Jin  (— ) 

®2 


®1  Jin  (^) 
®1 
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**1 

The  probe  calibration  factor,  — is  then  given  by, 

“2 


h 

No 


p a (r-> 


**2  ®2  In  (-^) 

^2 


(2.12) 


(v)  Curve  Fitting;  Modified  Prony's  Method 

The  experimental  data  obtained  has  magnitude  and  phase  so  that  it 
is  a complex  quantity.  Curve  fitting  to  the  experimental  data  is  done 
by  using  a modified  Prony's  method.  Prony's  method  [7]  fits  a set  of 
data  to  a series  of  exponentials.  By  assuming  a function  of  the  form 


f(z)  = A^e^^  + A2e~Y^ 


(2.13) 


because  of  the  TEM  mode  transmission,  only  two  terms  are  selected  from 
the  exponential  series  of  the  general  Prony's  method  with 

Y = (a  + j 3)  (2.14) 

where  a is  the  attenuation  constant,  3 is  the  propagation  constant. 

(2.15) 


f(z)  = 


where  A^^  and  A2  are  complex  constants,  A^^,  A2,  a,  3 are  to  be  determined. 
Let  Az  be  the  spacing  between  the  data  points  and 


^YAz 

(2.16) 

^-YASs 

(2.17) 

then 


Assume  and  U2  are  two  roots  of 


u + + Oq  = 0 


where  and  are  unknown. 


Equation  (2.18)  can  be  expanded  to 


fl  - Ai  + A^ 


2 2 
^3  ' V2 


c A k-1  ^ k-1 

^k  “ ■’■  ^2“2 


Since  the  above  set  of  equations  satisfy  equation  (2.19) 


^3  ^2“l  ^l“0  ° ° 


^4  ^3“l  + ^2“0  ° ° 
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(2.18) 


(2.19) 


(2.18) (a) 


(2.18) (b) 


(2.18) (c) 


(2.18)(d) 


The  unknowns  a^,  can  be  determined  from  the  above  equations.  Once 
Oq,  are  known,  the  roots  u^^,  U2  can  be  determined  from  equations 
(2.19).  Once  u^^  and  U2  are  known,  y he  obtained  from  equations 
(2.16)  or  (2.17)  or, 

(2.20) 


Y “ in(uj^)/Az 


-Y  - iln(u2)/Az 


J 


(2.21) 
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With  Y known  and  A2  are  then  determined  from  (2.18) (a)  and  (2.18)(b). 
Since  the  roots  and  U2  are  reciprocal  by  definition,  ^^2  ^ ‘ 

stitutlng  in  (2.19)  we  get 

“1^  + aj"*!  + “0  ■ “ 


~2  “iJT  “0  ■ 0 

“2  ^ 


(2.22) 


1 + aj^U2  + “0^2  “ ® 


(2.23) 


Comparing  (2.19)  and  (2.23)  we  obtain. 


“0  = ^ 


u + + 1 = 0 


(2.24) 


If  we  determine  then  the  roots  can  be  found  easily  in  the  following 


fl  + f2a^  + f3  = 0 


^2«1  “ -<^1  ^ ^3> 


^2  ^3“l  + ^4  = ° 


^3^  “ 


^3  ^4“l  + ^5  ' ° 


^4“l  “ “^^3  ^5^ 


With  three  datapolnts  one  can  easily  determine  from  which  and  U2 
can  be  determined. 

The  method  described  so  far  is  not  satisfactory  for  the  set  of 


data  where  experimental  noise  is  present.  A larger  number  of  datapolnts 
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Is  used  and  a least  square  curve  fitting  technique  is  employed. 

Let  MAz  be  the  spacing  between  datapolnts.  When  3M*N,  where  N 
is  the  total  number  of  datapolnts. 


f,  + f 0 ; f .,01  =-(f,  + fo  .,  ) 

1 nt+1  2m+l  ’ m+1  1 2m+l 


^2  ^nri-2“  ^2m+2  “ ° ’ ^m+2“  "^^2  ^2nri-2^ 


^3  ^iiH-3“  ^2nri-3  ° ’ ^mf3“  '^^3  ^2nri-3^ 


f + f,  a + f-  = 0 ; f,  a =-(f  + f_  ) 

m 2m  3m  2m  m 3m 


In  general 


[F]a  = -(G) 


= -[F^][G]  where  [F^]  Is  the  transpose  matrix  of  [F] 


Therefore  a = - [F^F]"^  [f’^]  [G] 


(2.25) 


The  residues  and  A2  can  be  determined  In  the  following  way 


f « A + A 

ti  + A^ 


, . yAz  , . -yAz 

^2  " V V 


fj^  - A^(e^^^)”  ^ 


[F]  - [G]  [A] 

[A]  » [g'^G]-^  [G^]  [F] 


(2.26) 
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With  A2,  Y]^  3nd  Y2  thus  determined,  a continuous  curve  may  then 
be  computed  from  equation  (2.15)  and  extrapolated  to  the  junction.  The 
value  of  and  in  equation  (2.12)  were  obtained  by  this  procedure. 
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CHARGE  PROBE. 


COAXIAL  SYSTEM  WITH  STEP  RADIUS  SYSTEM 


FREQUENCY 


DIAGRAM  OF  MEASUREMENT  SYSTEM 
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CHAPTER  III 


NUMERICAL  ANALYSIS 


(1)  Charge  Probe  as  a Scatterer;  Theory  and  Solution 


The  preceding  chapter  discussed  the  experimental  determination  of 
the  response  of  a charge  probe.  A theoretical  analysis  and  numerical 
solution  for  the  probe  response  and  calibration  factor  are  presented 
in  the  present  chapter.  For  the  purpose  of  this  analysis,  the  probe 
is  modeled  as  a monopole  scatterer  perpendicular  to  an  infinite  plane 
conductor  and  planar  image  theory  is  employed.  As  discussed  in 
Chapter  II  the  major  effect  of  the  cylindrical  conductor  upon  which 
the  probe  is  mounted  is  the  variation  of  the  E field  produced  by 
the  charge  on  the  test  structure.  In  the  treatment  of  the  probe  as  a 
scatterer  the  correct  dependence  of  the  incident  field  is  imposed. 

It  is  recognized,  however,  that  some  error  is  introduced  by  using 
planar  image  theory  rather  than  the  Green's  function  for  a source  ex- 
terior to  a cylindrical  conductor. 

From  the  specification  given  in  Chapter  II,  for  the  cable  from 
which  the  probe  is  constructed,  the  electrical  radius,  ka,  is  found  to 
be  .00036  for  an  operating  frequency  of  300  MHz.  It  is  therefore 


appropriate  to  make  the  usual  thin  wire  approximations. 
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(li)  Formulation  of  Straight  Wire  Integral  Equation 

The  charge  probe  of  length  slides  in  a charged  conductor  of 
radius  a^  in  a homogeneous  medium  characterised  by  (p,  e,  a = 0)  is 
shown  in  Fig.  (3.1).  The  z axis  of  a cylindrical  coordinate  system 
lies  at  the  center  of  the  charged  conductor  as  shown  in  Fig.  (3.1). 
The  probe  is  perpendicular  to  the  charged  conductor  and  extends  from 
p = to  p = For  convenience  a coordinate  transformation  is 

used  with  a new  variable  ^ = (p-a^^)  measured  along  the  probe  axis  as 
shown  in  Fig.  (3.2). 

The  current  and  charge  on  the  probe  produce  an  electric  field 

*8  S 

E (C)  on  the  probe  surface  which  has  a ^ component,  E^(?)  given  by. 


E®(Q  = ^ + k^)A_(C) 


(3.1) 


where  is  the  vector  potential  evaluated  on  the  probe  surface  and  k 
is  the  propagation  constant.  Since  A^(C)  is  produced  entirely  by  the 
current  on  the  probe,  it  is  given  by. 


i^(C’)K(c-c')d(;' 


C’=-H/2 


(3.2) 


where  I^(C)  is  Che  total  axial  current  on  the  probe  and 


“r  -Jk[(C-C’)^+  4a^  sin^(^)]^ 


K.,-?-)  = 2- 


^ r -jkt(^-^ 

J [(c-c’)^+ 

Ipo-TT 


4a^  sin^Y>r 
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where  a « radius  of  probe,  and  ii>  = angle  variable  of  a cylindrical 
coordinate  system  with  C as  the  axis. 

Applying  the  boundary  condition  that  the  tangential  electrical 
field  on  the  probe  surface  be  zero,  we  obtain, 

eJ(Q  + E®(g  = 0 

and  therefore 

eJ(C)  = -E®(C) 

Thus  the  Integro-dif f erentlal  equation  becomes. 


+ k^)  A^(C)  = E^(0 


(3.4) 
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local  behavior  at  C = 0 and  write. 

In  the  case  of  a scatterer, 

2 

^ ^ + k2)A^(;)  k 6(«  - eJ(C)  (3.5) 

As  explained  in  the  beginning  E^(^)  has  been  selected  to  be  propor- 
tional to  where  p = (C+a^).  Equation  (3.5)  with  (3.2)  is  the  final 
form  of  the  integral  equation  for  the  probe  current.  Equation  (3,5) 

has  been  solved  numerically  by  the  method  described  in  the  next  section. 

This  solution  technique  has  been  dealt  in  [8]. 

(iv)  Numerical  Solution  Technique 

In  this  section  a technique  for  solving  the  integro-dif f erential 

equation  is  presented.  This  technique  uses  piecewise  linear  testing 

and  pulse  expansion.  One  performs  the  testing  of  the  equation  to  be 

solved  with  elements  t^  of  the  selected  testing  set  before  the 

unknown  current  is  approximated  as  a linear  combination  of  the  elements 

1 of  the  basis  set C i 1. 
n \ n) 

Leading  to  the  numerical  solution  of  equation  (3.5),  one  equates 
the  corresponding  projections  of  the  two  sides  of  (3.5)  onto  the  space 
spanned  by  which  Implies  that  the  scalar  products  formed  by  t^ 

and  two  s''des  of  (3.5)  are  equated  for  each  m.  The  piecewise  linear 
testing  scheme  is  shown  in  Fig.  (3. A) 
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‘i?  ^ ’ *c  - v>  »•« 


(a) 


where  the  scalar  product  is 

^=y*  f 8*  (3.6)  (b) 

f ■)  C—H/2 

For  a real  testing  setVt  ^of  the  subdomain  type  considered  here  where 


+H/2 
f g*  dC 


t (C)  is  zero  outside 


Equation  (3.6)  becomes 


^ " <^m-r  ^m+l> 


Jii! 

k- 


ly  'jT  * " >*c«)  I?  +/  1^0)  2i„,,«(c)  t^(o  dt 


-^m+l 

C "C__n 


(C)  dc 


®“1  . 1 1 

I H c-  c„ 

Setting  t - A - ; 5-  - 

mm  A 

In  the  above,  the  equation  becomes 


(3.7) 

(3.8) 


m+l 


[l[l<  (C.)  * “7  *c  ■’'] 

'■'.-1 


^m+1 

/ 


* / "load  <C)  "S 


^m+1 

/ e; 


'm-1 


c-c. 


eJ(OA^  (C)  dC  (3.9) 


m-l 
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with  m - 0 + 1,  + 2,  + 3,  + 4,..,  + (M-1) 


(3. 10) (a) 


and 


“ • < W - Vi> 


(3.10) (b) 


For  a function  f sufficiently  smooth  over  the  Interval  (C  ,) 

m— 1 m+1 


one  employs  the  approximation 


^nri-l 


I 


f(C)  a"  (C)dC  - Af(C  ) 

m m 


(3.11) 


m-l 


The  final  equation  becomes 


<s>[ 


* * ^lo.d  «„> 


(3.12) 


where  5 is  the  kronecker  delta 
mn 


mn 


f 1,  m=n 

Y 0,  mj^n 


m = 0,  + 1,  ± 2,  + 3,  + 4 + (M-l) 


(v)  Expansion  of  Current  and  Formulation  of  Matrix 


The  current  is  represented  by  expanding  in  known  expansion 


functions  each  weighted  with  an  unknown  coefficient.  These  unknown 


coefficients  ultimately  determine  I^.  In  terms  of  a pulse  basis  set. 


^Pj^ , as  the  known  functions  and  the  complex  unknown  coefficients  I^, 
the  current  becomes 
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■ I „ 

n^-N 


(3.13) 


The;  pulse  arrangement  for  approximating  the  current  on  straight  wire 
Is  shown  In  Fig.  (3.5). 


The  vector  potential  can  be  written  as 


A.(C)  - I I A (C) 
^ n-N  " 


where 


H/2 

j 

C'=-H/2 


Pn(C')K(C-^')d^’ 


(3.14) 


(3.15) 


is  a partial  vector  potential  contributed  by  the  pulse  basis  element 
p^(C).  In  convenience,  the  following  quantities  are  defined: 


mn  , 2 
k 


^n+A/2  2 

45f  J - 2(1-  -P-  )K(e-6 


’n-A/2 


(3.16) (a) 


«=  A E"aj 
m 4 1° 

In  terms  of  the  above  definition  equation  (3.12)  becomes 


(3.16) (b) 


Z I = V 
mn  n m 


m=0,  +1,  + 2,.. + (M-1) 


(3.17) 
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The  boundary  condition  that  the  current  goes  to  zero  at  the 
ends  is  incorporated  by  directly  setting  P jj  “ Pjj  “ 0*  Equation 
(3.17)  can  be  conveniently  written  in  matrix  form  as, 

[Z][I]  = [V]  (3.18) 

Therefore  [l]  = [Z]~^[V] 

The  total  scheme  of  testing  and  expansion  is  shown  in  Fig.  (3.6). 

(vl)  Determination  of  Probe  Calibration  Factor 

Once  the  current  induced  in  the  load  can  be  found  as  a function 
of  a^^,  the  probe  calibration  factor  can  be  determined  as  follows. 

Let  be  the  current  Induced  in  the  load  when  the  probe  is  mounted 
on  the  larger  radius  tube  and  be  the  current  induced  in  the  load 
when  the  probe  is  mounted  on  the  smaller  radius  tube,  then  relating 
these  to  surface  charge  density 


ii  = 


(3.19) 


^2  ^2*^2 


(3.20) 


If  the  linear  density  of  charge  is  the  same  for  both  tubes  then 


qi  - q2 


or 


2lTaj^nj^=  2TTa2n2 


(3.21) 
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From  eq.  (3.19)  and  (3.20),  the  probe  calibration  factor 
obtained  as. 


Is 


(3.22) 


The  computed  value  of  this  ratio  for  different  probe  heights  for  two 
different  radii  of  cylinders  are  shown  in  Fig.  (4.19)  and  (4.20). 


Fig.  (3.2)  CHARGE  PROBE  AS  Fig.  (3.3)  LUMPED  LOAD  LOAD  AT  C-0. 

A SCATTERER. 


I 
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Fig.  (3.4)  TRIANGLES  IN  PIECEWISE  LINEAR  TESTING  SCHEME. 


CHAPTER  IV 


PRESENTATION  OF  DATA  AND  CONCLUSION 

(1)  Introduction 

Chapter  II  and  III  have  discussed  the  measurement  technique  and 
the  numerical  analysis  respectively.  It  Is  the  purpose  of  the 
present  chapter  to  present  the  data  resulting  from  both  approaches 
and  discuss  these  results. 

(11)  Measured  Data  for  the  Charge  Distribution 

As  explained  In  Chapter  II  on  measurement  the  charge  density 
measurement  was  made  for  3 charge  conditions  at  the  junction  near 
maximum  near  minimum  and  a condition  In  between  minimum  and  maximum. 

The  figures  (4.1)  to  (4.9)  present  the  charge  distribution  of 
the  step  radius  system  with  a^^*  5/16"  and  1/8"  for  3 different 
probe  heights  0.10",  0.08",  and  0.06".  For  each  probe  height, 
measured  data  for  3 different  charge  conditions  at  the  junction  are 
presented  along  with  the  curve  determined  from  the  modified  Prpny's 
method  for  the  magnitude  of  the  measured  data.  Similarly  the 
Figures  (4.10)  to  (4.18)  present  the  data  for  the  step  radius  system 
with  3/16"  and  1/8". 

(Ill)  Probe  Calibration  Factor 

The  probe  calibration  factor  (Nj^/N^)  as  a function  of  probe 
height  as  computed  from  the  measured  data  Is  shown  In  Fig.  (4.19)  for 
5/16"  and  aj”  1/8".  The  theoretical  probe  calibration 


factor 


obtained  by  the  technique  explained  In  Chapter  III  Is  also  shown  In 
Fig.  (4.19).  Similarly  Fig.  (4.20)  presents  the  probe  calibration 
factor  for  « 3/16"  and  a^  ■ 1/8". 

Examination  of  Figs.  (4.19)  and  (4.20)  reveals  differences  of  about 
8Z  In  the  experimental  values  of  N^/N^  for  different  charge  conditions 
are  the  junction.  With  the  assumption  that  the  probe  responds  to  the 
local  charge  density  at  Its  location  on  the  cylinder  should  be  a 

function  of  probe  height  and  cylinder  radii  only.  The  observed  effect 
of  the  position  of  the  junction  with  respect  to  the  standing  wave  pattern 
Is  believed  to  result  from  Inaccuracies  In  the  extrapolation  of  the  data 
to  the  junction.  The  extrapolation,  and  hence  determination  of  I^^/N2 
Is  accomplished  by  fitting  the  curve  to  the  measured  data.  Since  the 
data  Is  available  over  a distance  of  only  X/2  each  side  of  the  junction 
the  accuracy  In  determination  of  the  parameters  for  this  curve  Is 
limited.  For  these  reasons  the  data  obtained  for  the  case  of  maximum 
charge  at  the  junction  Is  the  most  reliable  since  the  slope  of  the  curve 
In  the  extrapolation  region  is  near  zero.  In  future  use  of  these  experi- 
mentally determined  probe  calibration  factors  greatest  weight  will,  there- 
fore, be  given  to  the  values  determined  for  the  maximum  charge  condition. 

Examination  of  Figs.  (4.19)  and  (4.20)  also  reveals  a consistent 

difference  between  the  measured  and  theoretical  values  of  N. /N^.  The 

1 2 

reasons  for  this  difference  are  not  fully  understood;  however.  It  must 
be  remembered  that  the  theoretical  solution  is  based  on  a planar  Image 
theory.  The  effects  of  the  cylindrical  geometry  are  the  most  likely 
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source  of  the  observed  differences  between  experimental  and  theoretical 
values.  Until  more  accurate  data  and  theoretical  solutions  are  available, 
the  experimental  value  of  Nj^/N2  will  be  used  as  the  most  reliable. 

(iv)  Conclusions 

The  results  of  this  Investigation  have  clearly  demonstrated  the 
effect  of  cylinder  radius  on  the  response  of  charge  probe  and  have  taken 
a considerable  step  toward  better  understanding  of  charge  measurements. 

The  probe  calibration  factors  obtained  can  now  be  used  to  correct  charge 
distribution  measurements  on  structures  containing  junctions  of  unequal 
radii  cylinders. 

Further  work,  however,  is  needed  for  complete  understanding  of  the 
response  of  monopole  probes  on  nonplanar  surfaces.  In  particular,  an 
effort  should  be  made  to  Include  the  cylindrical  image  plane  geometry 
into  the  theoretical  treatment  presented  here.  In  addition,  other  probe 

configurations  should  be  examined  and  evaluated  with  respect  to  their  | 

I 

response  on  nonplanar  surfaces.  | 
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Fig.  (4.1)  MAXIMUM  CHARGE  CONDITION  - PROBE  HEIGHT  = 0,10",  a,-  5/16 
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Fig.  (4.4)  MAXIMUM  CHARGE  CONDITION  - PROBE  HEIGHT  - 0.08",  a “ 5/16 
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Fig.  (4.5)  MEDIUM  CHARGE  CONDITION  - PROBE  HEIGHT  = 0.08",  a - 5/16' 
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MINIMUM  CHARGE  CONDITION  - PROBE  HEIGHT  = 0.08",  a = 5/16' 
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MAXIMUM  CHARGE  CONDITION  - PROBE  HEIGHT 
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Fig.  (4.8)  MEDIUM  CHARGE  CONDITION  - PROBE  HEIGHT  = 0.06",  a - 5/16 
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Fig.  (A. 11)  MEDIUM  CHARGE  CONDITION  - PROBE  HEIGHT  =0.10",  a - 3/16' 
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Fig.  (4.13)  MAXIMUM  CHARGE  CONDITION  - PROBE  HEIGHT  - 0.08",  a - 3/16' 
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MEDIUM  CHARGE  CONDITION  - PROBE  HEIGHT  - 0.08",  a - 3/16' 
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Fig.  (4.17)  MEDIUM  CHARGE  CONDITION  - PROBE  HEIGHT 
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MINIMUM  CHARGE  CONDITION  - PROBE  HEIGHT 


CALIBRATION  FACTOR  FOR  5/16"  to  1/8"  RATIO  OF  RADII 


Charge  Condition  at  the 


(A. 20)  PROBE  CALIBRATION  FACTOR  FOR  3/16"  to  1/8"  RATIO  OF  RADII 
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